Superconductivity in novel Ge-based skutterudites: {Sr, Ba}Pt4Gei2 



O 

o 

» 

o 



o 
o 
I 

Oh; 



I 

O 

o 



> 
oo 
m 

o 
d 

o 



X 



E. Bauer^, A. Grytsiv^, Xing-Qiu Chen^-*, N. Melnychenko-Koblyuk^, G. Hilscher^, H. 
Kaldarar\ H. Michor\ E. Royanian^ G. Giester^ M. Rotter^, R. Podloucky^ & P. Rogl^ 
^Institute of Solid State Physics, Vienna University of Technology, A-IO4O Wien, 



Institute of Physical Chemistry, University of Vienna, A-1090 Wien, 



Austria 
Austria 



Institute of Mineralogy and Crystallography, University of Vienna, Althanstrasse I4, A-1090 Wien, Austria and 
'^School of Materials and Metallurgy, Northeastern University, Shenyang, IIOOO4, P- R- China 

(Dated: February 2, 2008) 

Combining experiments and ab initio models we report on SrPt4Gei2 and BaPt4Gei2 as members 
of a novel class of superconducting skutterudites, where Sr or Ba atoms stabilize a framework 
entirely formed by Ge-atoms. Below Tc = 5.35 K, and 5.10 K for BaPt4Gei2 and SrPt4Gei2, 
respectively, electron-phonon coupled superconductivity emerges, ascribed to intrinsic features of 
the Pt-Ge framework, where Ge-p states dominate the electronic structure at the Fermi energy. 

PACS numbers: 74.70.Dd, 82.75.-z, 65.40.-b 



Cage-forming compounds such as zeolithes, fullerenes, 
clathrates or skutterudites have been proven not only 
of scientific but also of significant technological inter- 
est. The ability of these materials to accommodate guest 
filler species constitutes a wide range of varying chemi- 
cal and physical properties, comprising magnetic order, 
heavy fermion behavior, Fermi - and non-Fermi liquid 
features, or conventional or unconventional superconduc- 
tivity (SC). For recent reviews regarding the supercon- 
ducting properties of these classes of materials see Ref. 
[11, B S]. Most impressive, however, is the exceptional 
thermoelectric performance in some of the clathrates and 
skutterudites [J]. 

A structural-chemical description classifies skutteru- 
dites as cage compounds. The simple structure, however, 
acts as a prototype for a large class of compounds includ- 
ing binary as well as ternary and higher order represen- 
tatives. Hitherto, cage forming elements of skutterudites 
are essentially based on volatile and/or toxic pnictogens 
(X = P, As, Sb). Building blocks in the framework are 
8 tilted octahedra per unit cell enclosing two icosahedral 
cages. Each of the octahedra is centered by a transition 
metal atom like Co or its homologues (Rh, Ir). 

Motivated by the known manifold of interesting phys- 
ical properties among skutterudites, we searched for 
ternary and/or higher order compounds exploiting the 
combination of (i) a high density of d-states of a platinum 
group metal at the Fermi level with (ii) a rigid frame- 
work of rather covalently four-bonded atoms such as Si 
and/or Ge. Thereby a novel family of Ge based skut- 
terudites, SrPt4Gei2 and BaPt4Gei2 has been identified. 
The aim of the present work is a study of the stability and 
the characterization of bulk properties by means of resis- 
tivity, magnetization, specific heat and band structure 
calculations. For both compounds, the low temperature 
behavior is dominated by the appearance of supercon- 
ductivity around 5 K. 

{Sr, Ba}Pt4Gei2 were prepared from stoichiometric 
amounts of high purity elements by argon arc melting 
and subsequently heat treated in evacuated quartz cap- 
sules at 800°C for two weeks. Phase purity and lattice 



parameters were checked by EMPA and x-ray diffrac- 
tion. Bulk properties of these novel skutterudites were 
obtained by a number of standard techniques, details are 
summarized in Ref 5*1 . Density functional theory (DFT) 
was applied using the Vienna ab initio simulation package 
(VASP) [aja with a fully relativistic spin-orbit coupling 
approach [oj. The Brillouin zone was sampled with 
5x5x5 Monkhorst-Pack /c-point grids. The exchange- 
correlation functional was treated within the local den- 
sity approximation. 

The crystal structure of {Sr, Ba}Pt4Gei2 was deter- 
mined from Kappa-CCD single crystal X-ray data and 
found to be cubic, isotypic with the filled skutterudite 
type LaFe4Sbi2 [1C| . Structure and lattice parameters 
are collected in Tabled] Occupation factors were refined, 
corresponding to a full occupancy of the Pt and Ge sub- 
lattices. Although not revealed from single crystal re- 
finement, minor deviations (up to 3 %) from full occu- 
pancy are possible for the Ba and Sr atoms. Since the 
size of the Ge- framework is significantly smaller than the 
corresponding Sb-framework, effective bonding between 
the framework cages (two icosahedra per unit cell) and 
the Ba(Sr)-center atoms is ensured. As a consequence 
of this stronger bond between cage and guest atom, we 
observe very regular thermal atom displacement factors 
(ADP) on all atoms. The temperature dependencies of 
ADP's in the temperature region from 100 to 300 K re- 
veal for all atoms similar trends: typical rattling modes 
caused by the heavy guest atoms in Sb-based skutteru- 
dites are absent in {Ba, Sr}Pt4Gei2. The structural pa- 
rameters as derived from DFT calculations are in ex- 
cellent agreement with the experimental data. In order 
to analyze the thermodynamical stability of XPt4Gei2 
(X=Ba,Sr) compounds we also calculated the total en- 
ergy for a hypothetical compound Pt4Gei2. The bonding 
energy Ex for guest atom X is obtained from the relation 
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notes the corresponding total energy of the compound or 
elemental solid in its equilibrium state. The values i?Ba 
= -3.24 eV and i?sr = -3.38 eV emphasize the stabilizing 
effect of the Ba and Sr guest atoms. 
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FIG. 1: (Color online) Temperature dependent electrical re- 
sistivity p of BaPt4Gei2 (a) and SrPt4Gei2 (b). Both insets 
show the field dependence of p(r). 



Superconductivity is found for BaPt4Gei2 and 
SrPt4Gei2 from resistivity measurements using polycrys- 
talline samples at Tc w 5.3 K and 5.1 K, respectively 
(compare Fig. [Ha,b)). The normal state regions, T > Tc 
of these ternary compounds do not behave like simple 
metals, since the standard model of the electrical resistiv- 
ity of metallic systems, i.e., the Bloch-Griineisen formula 
is not applicable. Such observations were made in many 
superconducting materials and may be attributed to a 
substantial electron-phonon interaction strength, respon- 
sible for the formation of Cooper pairs in conventional 
superconductors. Rather, the overall p(r)-dependence of 
both skutterudites follows the model of Woodward and 
Cody which initially was applied to A15 supercon- 
ductors such as NbaSn. Least squares fits to this model 
are shown as solid lines in both figures, revealing reason- 
able agreement for characteristic temperatures Tq = 123 
and 121 K, respectively. The differences of the residual 
resistivities may correspond to small differences in the 
filling of the 2(a)-sites by Ba or Sr. In case of the smaller 
atom Sr, filling seems to be more complete. The applica- 
tion of a magnetic field suppresses superconductivity of 
BaPt4Gei2 at fields above 2 T, while for SrPt4Gei2 the 
upper critical field -ffc2(0) ~ 1 T (insets of Fig. [1]). 

The susceptibility x exhibits a rather sharp transition 
at r = 5.3 and 5.1 K [inset. Fig. [HJb)], dropping from an 
initially small positive value to the diamagnetic value of 
[— l/(47r)] for zero field cooling, which corresponds to a 
full Meissner Ochsenfeld effect. Magnetization measure- 
ments performed at various temperatures and magnetic 
fields up to 6 T evidence type II superconductivity. Up- 
per critical field values are summarized in Fig. [31 

The heat capacity, Cj,, of {Sr, Ba}Pt4Gei2 is plotted 
in Fig. m as Cp/T vs. T for zero and 3 T. For sake 
of clarity, measurements taken at intermediate field val- 
ues are not shown here. The jump of Cp(T) below 6 K 
evidences bulk superconductivity in both cases. An ide- 
alization of the heat capacity anomaly under the con- 
straint of entropy balance between the superconducting 
and the normal state yields Tc = 5.35 and 5.1 K for 



TABLE I: Normal state and superconducting properties of 
BaPt4Gei2 and SrPt4Gei2 which crystallize in the cubic skut- 
terudite structure: space group Im3, (No. 204); Ba and Sr 
are at the 2(a) (0, 0, 0) sites, Pt at the 8(cJ (i, i, \) sites, 
and Ge at the 24^ (0, y, z) sites. Ueq is a mean value of the 
atomic displacement ellipsoid. 



property 


IDar^t4ljrCl2 
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lattice parameter a @300 K [nm] 


0.86928(3) 


0.86601(3) 


Ge 24^ site: y 


0.15302 


0.15197 
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Rf2 = Y:\Fo -F^\/Y.F^ 
f/e,({Ba,Sr}) [nm^] 
f/e,(Pt) [nm^] 


0.019 


0.018 


0.000066(2) 


0.000118(3) 


0.000062(1) 


0.000069(1) 


f/eg(Pt) [nm^] 


0.000092(2) 


0.000097(2) 


transition temperature Tc [K] 


5.35 


5.10 


upper critical field jj,oHc2 [T] 


1.8 


1 


thermod. critical field poHc [mT] 


53 


52 


Fermi velocity vf [m/s] 


52500 


67000 


coherence length ^ [nm] 


14(1) 


18(1) 


penetration depth A [nm] 


320(10) 


250(10) 


G.L. parameter k 


24(1) 


14(1) 



the Ba and Sr based compounds, respectively. Assuming 
that the specific heat of metallic compounds at low tem- 
perature follows Cp = 7T -I- f3T^, (7 is the Sommerfeld 
coefficient and P is proportional to the Debye tempera- 
ture 6u), least squares fits reveal 7 = 42 mJ/molK^ and 
9d = 247 K (Ba) and 7 = 41 mJ/molK^ and 9d = 220 K 
(Sr). It is worth to be noted that the Debye tempera- 
ture of BaPt4Gei2 is larger than that of SrPt4Gei2. In 
general, however, materials with smaller masses exhibit 
larger Debye temperatures. This anomaly observed may 
correspond to the fact that, while the volume of the unit 
cells of both compounds differ by only 1%, the atomic 
volumes of Sr and Ba differ by about 12 %. This causes 
a weaker bonding of Sr to the framework, hence a weaker 
force constant may result in lower values of Od- 
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FIG. 2: (Color online) (a) Temperature dependent specific 
heat Cp of SrPt4Gei2 and BaPt4Gei2 plotted as Cp/T vs. T 
for and 3 T. (b) SemHogarithmic plot CeshTc vs. Tc/T. 
The inset shows the low temperature and low field suscepti- 
bility M/H for both compounds. 
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FIG. 3: (Color online) (a) Temperature dependent upper crit- 
ical field ^oHc2 of BaPt4Gei2 and SrPt4Gei2 as obtained 
from resistivity, magnetization and specific heat measure- 
ments. The dashed and the dashed-dotted lines represent 
the WHH model. 



Taking into account the McMillan model ^] allows 
calculation of the dimensionless electron-phonon cou- 
pling constant A. Setting the repulsive screened Coulomb 
term /u* k, 0.13, yields in both cases A « 0.7. This refers 
to superconductors well beyond the weak coupling limit. 
In comparison, /i* of different cage forming compounds 
covers a range from « 0.1 to « 0.3 [H, d, [1]. 

The jump of the specific heat ACp/T(T = T^) w 
58mJ/molK^ (Ba) and ~ 57mJ/molK^ (Sr) allows calcu- 
lation of ACp/(7„Tc) ~ 1.35, which is near to the figure 
expected from BCS theory [ACp/(7Tc) w 1.43]. As the 
magnetic field strength increases (not shown here), both 
the transition temperature and the anomaly right at Tc 
are suppressed, constituting the phase diagram in Fig. [31 

The superconducting gap A(0) can be derived from a 
comparison of the modified BCS expression, Ces{T) — 
8.57rcexp(-0.82A(0)/A:BT) with the experimental data 
depicted in a semi- logarithmic plot Ces/lT,. vs. T^/T 
in Fig. [Ifb) where BaPt4Gei2 and SrPt4Gei2 exhibit 
for Tc/T > 2 an exponential temperature dependence 
indicating a ratio A(0)/A;Brc ~ 1.8 in close agreement 
with the BCS value Abcs(O) = l.lQk^Tc. 

The thermodynamic critical field is calculated from the 
free energy difference between the superconducting and 
normal state: AF(T) = F„-F^ = fioH^{T)/2, where F„ 
and Fg are evaluated from the specific heat in the nor- 
mal and superconducting state, respectively. The values 
obtained are /lo-ffc(O) ~ 53 (2) and 52 (2) mT for Ba and 
the Sr-based compound, respectively. 

Fig. [3] displays the temperature dependent upper crit- 
ical field fJ.oHc2 of {Sr, Ba}Pt4Gei2 as deduced from field 
dependent resistivity, magnetization and heat capacity 
measurements. The slopes of the upper critical field 
d{fioHc2)/dT = fJ.oH'i.2 are collected in Table U yielding 
slightly larger values deduced from magnetization and re- 



sistivity data than those from specific heat, which may 
be attributed to pinning and surface effects. /io-ffc2 of 
BaPt4Gei2 is larger than /ioi?c2 of SrPt4Gei2. 

Primarily, two mechanisms are responsible for the lim- 
ited value of fJ.oHc2- orbital pair breaking and Pauli 
limiting. Werthamer et al. |lj|| derived an expression 
( WHH model) for the upper critical field ^oHc2 in terms 
of orbital pair-breaking, including the effect of Pauli spin 
paramagnetism and spin-orbit scattering. A comparison 
of the experimental results with the WHH model is based 
on two parameters, a, the Pauli paramagnetic limitation 
{Maki parameter) and A^o describing spin-orbit scatter- 
ing. If the atomic numbers of the elements constituting 
the material increase, Xso is expected to increase as well. 

The Maki parameter a can be estimated via 
the Sommerfeld value 7 and po [Tsj . i.e., a — 
(3e^?i7Po)/(2TO7r^fc^) with e the electron charge and m 
the electron mass. Taking the experimental po and 
7 yields a = 0.18 for BaPt4Gei2 and a = 0.14 for 
SrPt4Gei2. We have adjusted the WHH model to the ex- 
perimental data (dashed and dashed-dotted lines in Fig. 
[3l), revealing A^o ~ 15 for all data-sets. 

The thermodynamic and the upper critical field are 
used to calculate the ratio of the penetration depth 
Agl(O) to the coherence length ^gl(O) via Abrikosov's re- 
lation AGL(0)/fGL(0) = kgl{0) = H,2/[V2H,{Q)] yield- 
ing the Ginzburg-Landau parameter kql ~ 24 (2) and 
14 (2) for the Ba and Sr-based compounds. The absolute 
values of the coherence length ^0 and the penetration 
depth A(0) can be evaluated via the isotropic Ginzburg- 
Landau-Abrikosov-Gor'kov (GLAG) theory. Values ob- 
tained are presented in Table HI 

The Fermi velocity vp can be calculated from the 
Fermi surface area Sg as shown by Orlando et al. p^ . 
Values are given in TablelH Combining Ss and po, a mean 
free path kr of about « 1.0 • 10"® m and w 1.4 • 10~® m 
for the Ba and Sr compound, respectively, is derived. 

The fact that /io-f^c2 of BaPt4Gei2 is larger than that of 
SrPt4Gei2 can be understood in terms of the Ginzburg- 
Landau theory. Hake jl^ and Orlando et al. [l3| derived 
a model equation for poH'^2 which primarily depends on 
two parameters: on the inverse of the square of and 
on the inverse of Itr- Taking into account the parameters 
deduced from the above analyses explains in both cases, 
at least qualitatively, differences of /ioi?c2- While the 
decrease of fi? from the Sr to the Ba case can be conceived 
by an increase of the unit cell volume {vp o(: (N/V)^^^, 
for free elelctrons), the decrease of the mean free path 
corresponds to the increase of the residual resistivity from 
SrPt4Gei2 to BaPt4Gei2. 

From kr/S, ~ 1 we classify {Ba, Sr}Pt4Gei2 as a dirty 
limit superconductor, and k of the order of 10 to 20 refers 
to a pronounced type II superconducting behavior. 

The DOS in Fig. 4 reveal the individual contributions 
of X=(Ba,Sr), Pt and Ge atoms showing that Ge states, 
which are of p— like character, are dominating at Fermi 
energy. The Ge states hybridize with Pt 5c?-like states by 
which the spin-orbit coupling effect is transmitted to the 
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FIG. 4: (Color online) DOS for BaPt4Gei2 (a,b) and 
SrPt4Gei2 (c,d) as derived from DFT calculations. Fermi en- 
ergy is situated at zero energy. The upper half (a),(c) shows 
the density of states for a fully relativistic calculation, which 
includes spin-orbit (SO) coupling. The lower half (b),(d) of 
the figure presents the DOS of a standard non-relativistic cal- 
culation. The total DOS is decomposed into the local DOS 
representative for one Ba(Sr), four Pt, and twelve Ge atoms. 
The relativistic effect is small but important, because it shifts 
the Fermi level of the non-relativistic DOS almost into the 
maximum of the relativistic DOS. 



DOS at E^. Consequently, by making use of the electron 
localization function technique [? ], expressive Ge-Ge co- 
valent bonds and typical metallic Pt bonding are being 
deduced. The metallic features of the DOS around Ep 
convincingly confirm that the Zintl concept no longer ap- 
plies to {Sr, Ba}Pt4Gei2, while its applicability is quite 
obvious among pnictogen-based skutterudites. 

The total DOS at Ep can be compared with the Som- 



merfeld value of the specific heat 7 = ^Tr'^N{Ep)kg. The 
experimental values slightly larger than 40 mJ/molK^ 
for the Ba and Sr compounds, along with an electron- 
phonon enhancement factor X^p — 0.7 (estimated via 
the McMillan formula [13]) would require a bare DOS 
equivalent to « 25 mJ/molK^, which compares favorably 
with the DOS calculations involving spin-orbit coupling: 
31 mJ/molK2 for the Ba {N{Ep) = 13.2 states eV^^ 
f.u^i) and 28 mJ/molK^ for the Sr compound iN{Ep) 
= 12.1 states eV'^ f.u-^). Without spin-orbit coupling, 
the specific heats, 7, are 19.5 and 21.0 mJ/molK^ for Ba 
and Sr, respectively. This implies the importance of the 
relativistic effects. 

In summary, superconducting {Sr, Ba}Pt4Gei2 are the 
first skutterudites where the framework in the structure 
is entirely built by Ge-atoms. DFT calculations proved 
that X=(Ba,Sr) guest atoms strongly stabilize the com- 
pounds. Most strikingly, the calculated DOS around 
Fjp is composed of hybridized Ge 4p-like and Pt 5c?-like 
states, and it has a sharp peak with its maximum very 
close to Ep. The influence of the guest atoms (Ba or 
Sr) on superconductivity, however, may be ruled out due 
to the fact that (i) the Ba- or Sr-likc DOS around Ep 
is negligible and (ii) the DOS around Ep for the hypo- 
thetical X-free Pt4Gei2 framework is very similar to the 
one of XPt4Gei2 (X = Ba and Sr). Hence, superconduc- 
tivity appears to be an intrinsic property of the Pt-Ge 
cage-forming structure. This conclusion is in line with 
the slightly smaller value of Tc observed in SrPt4Gei2, 
in marked contrast to the isotope effect, where lighter 
masses would originate larger SC transition tempera- 
tures. 
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